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Ab initio lattice dynamics simulations and inelastic neutron scattering spectra for 
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We have performed extensive ab initio calculations to investigate phonon dynamics and their 
possible role in superconductivity in BaFe2As2 and related systems. The calculations are compared 
to inelastic neutron scattering data that ofTer improved resolution over published data [Mittal et al., 
PRB 78 104514 (2008)], in particular at low frequencies. Effects of structural phase transition and 
full/partial structural relaxation, with and without magnetic ordering, on the calculated vibrational 
density of states are reported. Phonons are best reproduced using either the relaxed magnetic 
structures or the experimental cell. Several phonon branches are affected by the subtle structural 
changes associated with the transition from the tetragonal to the orthorhombic phase. Effects of 
phonon induced distortions on the electronic and spin structure have been investigated, ft is found 
that for some vibrational modes, there is a significant change of the electronic distribution and 
spin populations around the Fermi level. A peak at 20 meV in the experimental data falls into the 
pseudo-gap region of the calculation. This was also the case reported in our recent work combined 
with an empirical parametric calculation [Mittal et al., PRB 78 104514 (2008)]. The combined 
evidence for the coupling of electronic and spin degrees of freedom with phonons is relevant to the 
current interest in superconductivity in BaFe2As2 and related systems. 

PACS numbers: 74.25.Kc,74.70Dd,78.20.Bh,78.70.Nx 
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I. INTRODUCTION 



Every time a new class of superconducting material 
is discovered another door is opened to an improved 
understanding of the phenomenon of high temperature 
superconductivity. This holds, in particular, when 
transition temperatures beyond some 20 K are reached. 
Examples are the research activities triggered by the 
discovery of superconductivity in doped fuUerenes, 
boron-carbides and MgB2. Despite the undeniable 
new insights procured by these studies, the question 
of high-temperature superconductivity continues to 
elude description by physical models. Oxypnictides, a 
family of superconductors with rather high transition 
temperatures that share structural and electronic two- 
dimensonality with copper-oxides, offer a new class of 
materials to study high temperature, superconductivity. 
Superconductivity has been evidenced in iron- 
oxypnictide LaOFeP— with Tc ~ 4 K, increasing to 
~ 7 K with F doping^. A similar transition temperature 
has been measured for nickel-oxypnicitde LaONiP'^ 
(Tc - 4 K). An unexpectedly high-T^ (T^ - 26 K) 
was recently discovered in the F-doped iron-oxypnicitde 
La[0i_2;Fa;]FeAs^i^. The parent compound LaFeAsO 
offers a unique picture of the chemical bonding; dom- 
inantly ionic in the LaO layers and covalent in the 
FeAs layers. Thus charge transfer can be assumed to 
occur according to (LaO)"'' (FeAs) ^. It is established 
that the superconductivity in LaFeAsO originates 
mainly from structural and electronic properties of 



the FeAs layers. Similar structural and electronic 
features have been observed in the ternary, oxygen-free, 
iron-arsenide BaFe2As2. The only difference between 
the two compounds is that the FeAs layers, which are 
in principle electronically identical in both cases, are 
separated in BaFe2As2 by Ba atoms instead of LaO 
layers. In this case the charge transfer occurs according 
to BaQ~^(FeAs)~. BaFe2As2 becomes superconducting 
under pressure at 29 K® and exhibits a phase tran- 
sition around 140 K involving a structural transition 
(distortion) from tetragonal (high temperature) to 
orthorhombic (low temperature)^^. BaFe2As2, featuring 
the simplest possible buffer in between the FeAs-layers, 
is therefore ideal for studying the coupling of nuclear 
and electronic degrees of freedom. 

Phonons and electron-phonon coupling have recently 
been studied computationally in the parent compound 
LaOFeAsM. It is generally concluded that conventional 
electron-phonon coupling is insufficient to explain 
superconductivity in the family of Fe-As-based super- 
conductors, while some reports in literature indicate 
the possibility of non-conventional electron-phonon 
coupling. The example of MgB2 showed that phonon 
coupling in conjunction with particular Fermi surface 
nesting effects is a viable way to describe transition 
temperatures of up to 40 K— . If phonons do not 
directly mediate the electron pairing interaction, other 
excitations are required and spin fluctuations have been 
considered and measured in BaFe2As2ii. 
Experimentally, phonons in Fe-As superconductors have 
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FIG. 1: (Color online) Top panels: dispersion relations and vibrational DOS for tetragonal (SI) and orthorhombic (S2) phases. 
For the tetragonal phase the high symmetry points are N=(l/2,0,0), X=(0, 0,1/2) and Z=(l/2, 1/2,-1/2). For the orthorhombic 
phase they are L=(l/2,0,0), Y=(l/2,0,l/2) and Z=(l/2,l/2,0). Middle panels: Weighted, partial, vibrational DOS and the 
generalised vibrational DOS for the two phases. Experimental, generalised vibrational DOS, measured at two temperatures. 



been investigated by inelastic neutron scattering (INS) 
in fluorine doped LaOFeAs^^ and BaFe2As2^. Thermal 
neutrons were used, allowing spectra to be measured in 
neutron down-scattering from 4K up to room tempera- 
ture, but this spectrometer configuration tends to limit 
energy resolution for low frequency modes, compared to 
using cold-neutrons in an up-scattering configuration. 
In LaOFeAs, experimental data was compared with 
published, ab initio phonon calculationsH. The agree- 
ment between experiment and calculation is reasonable, 
corresponding band frequencies being within a few 
meV, but this level of inaccuracy leads to strong peaks 
at 30 and 40 meV in the experimental spectra being 
compared to a single peak at 35 meV in the calculation. 
In the case of BaFe2As2 the experimental resolution was 
not sufficient to see evidence of the structural phase 
transition at 140 K in the phonon spectrum. In order to 
better evaluate coupling between electronic and nuclear 
structure and phonons, higher resolution INS data is 
required that can be rigorously compared with ab initio 



calculations. 

In this context, we present a combined theoretical and 
experimental investigation of the lattice dynamics of 
BaFe2As2, including coupHng to electronic and spin 
degrees of freedom. We report higher resolution INS 
spectra, above and below the phase transition tem- 
perature, which complement our published data. The 
generalised vibrational density of states is calculated 
for both phases using ab initio methods and allows the 
atomic contributions to the phonons to be evaluated. 
Phonon calculations have been repeated on experi- 
mentally determined crystal structures and on relaxed 
geometries, including spin polarisation, in order to 
identify critical, structural degrees of freedom. Finally 
the effect of phonons on the spin-resolved, electronic 
density of states is investigated. 



TABLE I: Labels for the different ab initio simulations presently reported. 



Label 



Structural Data 



Phase 



SI 
S2 
S31 
S32 

S41 
S42 



experimental 

experimental 

fully relaxed 

partially optimized (only z(As) is optimized) 

fully relaxed + magnetic ordering 



tetragonal 
orthorhombic 
orthorhombic 
orthorhombic 
orthorhombic 



partially optimized (only z(As) is optimized) + magnetic oredering orthorhombic 



II. EXPERIMENTAL DETAILS 

The neutron spectroscopy measurements were per- 
formed on the cold time-of-flight spectrometer IN6 (Insti- 
tut Laue Langcvin, Grenoble, France), in up-scattering 
mode using an incident wavelength of 5.1 A and time- 
focussing at 4 meV giving very good resolution below 20 
meV. Compared to our previous worki^, where the ther- 
mal neutron time-of-flight spectrometer IN4 was used in 
down-scattering, we have better resolution but cannot 
access the full spectrum below a temperature of 140 K. 
The measurements from the two instruments are in this 
respect fully complementary and both structural phases 
can just be measured on IN6. Two grams of sample of 
polycrystalline BaFe2As2 were prepared and character- 
ized as reported elsewherel^. The sample was held in 
a thin Al-foil container and thermalized within a nitro- 
gen cryo-loop for measurements at 140 and 300 K. The 
container scattering was subtracted and the detectors 
calibrated with a vanadium standard. After correction 
for energy dependent detector efficiency, the one-phonon, 
generalized phonon DOS (GDOS) was obtained from the 
angle-integrated data (17 to 114°) using the incoherent 
approxinration-^'. Multi-phonon terms to the scattering 
were determined self-consistently in an iterative proce- 
dure. Figure 1 (bottom panel) shows the measured spec- 
tra. The temperature of the structural transition is re- 
ported to be 140 K— or 142 K— , with peak intensities 
starting to change significantly at temperatures as high 
as 150 K. The 140 K spectrum is therefore measured very 
close to the transition temperature but it shows some no- 
table differences with the 300 K spectrum, like the hard- 
ening of the highest frequency vibrations and the band 
at 15 meV. These trends tend to be reproduced by the 
calculations (see below). Bigger differences between the 
spectra of the two structural phases are not expected 
since the data measured over a wide temperature range, 
at lower resolution, on IN4 did not reveal any marked 
changesi^. 



mented in the Vienna ab initio simulation package 
(VASP)2ii22. The GGA was formulated by the Perdew- 
Burke-Ernzerhof (PBE)^-2j^ density functional. The 
Gaussian broadening technique was adopted and all re- 
sults are well converged with respect to fc-mesh and en- 
ergy cutoff for the plane wave expansion. 
Experimentally refined crystallographic data in both the 
low-temperature (orthorhombic) and high-temperature 
(tetragonal) phases were taken from referencel-. These 
structures were used to calculate the GDOS and disper- 
sion relations for both the tetragonal and orthorhombic 
phases. Hereafter, the corresponding simulations are la- 
belled SI and S2, respectively. The orthorhombic phase 
was then taken for a case study for which full (S31) and 
partial (S32) structural relaxation were performed. By 
partial relaxation we mean only the As z position is opti- 
mized. The effect of the observed magnetic ordering— 
in the orthorhombic phase was taken into account in 
further full (S41) and partial (S42) optimizations. Ta- 
ble I provides a summary of these simulations. In the 
magnetic structure the Fe moments are parallel to the 
longer a-axis. Spins are aligned ferromagnetically along 
the shorter b-axis and antiferronragnetically along the a 
and c axeaii. 

In the lattice dynamics calculations, in order to deter- 
mine all inter-atomic force constants, the supercell ap- 
proach has been adopted^. Therefore, the single cells 
were used to construct (3*a, 3*b, c) and (2*a, 2*b, c) 
supercells for the tetragonal and orthorhombic cases, re- 
spectively, a and b being the shorter cell axes. The for- 
mer contains 18 formula-units (90 atoms) and the lat- 
ter contains 16 formula- units (80 atoms). Total energies 
and inter-atomic forces were calculated for the 12 and 
18 structures (in the tetragonal and orthorhombic su- 
percells, respectively) resulting from individual displace- 
ments of the three symmetry inequivalent atoms, along 
the three inequivalent Cartesian directions (±x, ±y and 
±z). The 15 phonons branches corresponding to the 5 
atoms in the primitive cell, were extracted in subsequent 
calculations using the Phonon software^^. 



III. COMPUTATIONAL DETAILS 



Electronic and ionic first principle calculations were 
performed using the projector-augmented wave (PAW) 
formalismi^ of the Kohn-Sham DFTASi^ at the gen- 
eralized gradient approximation level (GGA), imple- 



TABLE II; Experimental and optimized lattice parameters and fractional As z-coordinates for both structural phases. NM, F 
and Mobs label non-magnetic, ferromagnetic and observed magnetic ordering. Values between brackets indicate optimized As 
z position obtained from partial structural relaxations. 







Orthorhombic 


phase 




Tetragonal phase 




Expia 


Mobs 




NM/F 


Expia 


Mobs 


NM/F 


a 


5.6146 


5.6611 




5.6026 


3.9625 


3.9746 


3.9553 


b 


5.5742 


5.5940 




5.5998 


3.9625 


3.9746 


3.9553 


c 


12.9453 


12.8793 




12.6302 


13.0168 


12.7154 


12.6849 


z(As) 


0.3538 


0.3512 (0.3510) 


0.3456 (0.3446) 


0.3545 


0.3469 (0.3467) 


0.3455 (0.3442) 



IV. RESULTS 



Generalized phonon density of states and 
dispersion relations 



Figure 1 shows the calculated dispersion curves and the 
GDOS for both the tetragonal (high-temperature: SI) 
and orthorhombic (low-temperature: S2) phases using 
the experimental structures"'^^. In order to compare with 
experimental data, the calculated GDOS was determined 
as the sum of the partial DOS weighted by the atomic 
scattering cross sections and masses: GDOS=^j jf- 
pDOSj ijf- = 0.21 (Fe), 0.07 (As), 0.02 (Ba); i={Fe, As, 
Ba}). Both the full (GDOS) and partial gdos (pDOS) 
are included in figure 1. The comparison of the calculated 
and measured spectra is generally good. The acoustic re- 
gion extends in both phases to about 9 meV. The lowest 
lying optic modes give rise to intensity starting at about 
12 meV. As can be seen from the phonon relations, the 
corresponding dispersion sheets are not flat. While at 
the F-point these low frequency, optic modes correspond 
to rigid displacements of the FeAs layers with respect to 
the Ba spacers, distortions within the FeAs layers are 
present in the zone boundary modes. F-point modes 
above 18 meV distort the FeAs layers and the Raman 
active ones are illustrated in Figure 2(c). The calculated 
vibrational spectrum is clearly separated into three bands 
via two pseudo-gaps, which arise from the strongly dis- 
persive sheets in the regions around 20 and 28 meV. 
The calculated phonon dispersion is sensitive to struc- 
tural details, particularly for the optic bands in the inter- 
mediate frequency range from 10 - 20 meV. Flat bands 
lead to a clear two-peak structure in the GDOS at 13 
and 18 meV in the orthorhombic phase. These peaks are 
found at 13 and 16 meV in the tetragonal phase in good 
agreement with peaks in the corresponding experimental 
data. In particular, the experimentally observed peak at 
16 meV in the orthorhombic phase softens to 15 meV in 
the tetragonal phase. Otherwise at the highest frequen- 
cies, bands at 24 and 29 meV in the tetragonal phase 
harden to 25 and 31 meV in the orthorhombic phase, fol- 
lowing the trends observed experimentally. It is the new, 
higher resolution INS data that allows the validity of the 
phonon calculations, with respect to the structural phase 
transition, to be established. 
However there is a peak observed at 20-21 meV which is 



absent in the calculated spectra for both phases. This 
constitutes the only significant discrepancy between ex- 
periment and calculation. As can be seen in the partial 
DOS, the dispersion branches in the region under ques- 
tion represent motions of Fe and As. 



B. Effect of structural relaxation 

A generally good description of the phonon spectra has 
been obtained by using the experimental structures for 
the ab initio calculations. However, experimental struc- 
tures are not exactly the equilibrium structures that are 
typically used for phonon calculations. Phonon calcula- 
tions have therefore been repeated in the orthorhombic 
phase for which full (S31 - atomic z co-ordinate of As 
atom plus lattice constants) and partial (S32 - atomic 
co-ordinate only) structural relaxation were performed. 
In Figure 2(a) the corresponding spectra are compared 
with the calculation performed on the experimental cell 
(S2) and the experimental spectrum. The S31 and S32 
spectra are almost identical meaning that the lattice con- 
stants alone have only a rather minor influence on the 
phonon spectra. The lattice constant which changes most 
on full relaxation is the c parameter, which determines 
the spacing between Ba and Fe-As layers (see Table II). 
The highest frequency modes, above 24 meV react most 
significantly to these structural changes. For example, 
the band at 25 mcV splits in both S31 and S32, whereas 
a much weaker splitting is observed in the experimental 
data. The agreement with the experimental data is worse 
than for the calculation performed on the experimental 
structure, due to the hardening of the highest frequency 
vibrations. The only internal structural parameter is the 
z coordinate of the As atom, which deviates by 0.1 A 
from the experimental value upon relaxation. This is 
cause of the disagreement at high frequency. 



C. Effect of magnetic ordering 

Further full (S41) and partial (S42) structural relax- 
ation in the orthorhombic phase were performed taking 
into account the observed magnetic orderingii (see Fig- 
ure 2(b)). This leads to a description of the GDOS which 
is significantly better than that shown in Figure 2(a) for 
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FIG. 2: (Color online) Measured and simulated GDOS in the orthorhombic phase (a) for relaxed non-magnetic structures (S31 
and S32) and (b) for relaxed structures including the known magnetic ordering (S41 and S42). In both cases, the calculation 
based on the experimental, orthorhombic structure S2 is shown for comparison. Labels B, C and E in (b) indicate specific 
phonon modes which will be considered to discusse their effect on the electronic structure of BaFe2As2(cf Figure 3). The full 
displacement patterns of the six Raman active modes (A: 18 meV, B: 18 meV, C: 26 meV, D: 29 meV, E: 33 meV and F: 34 
meV ), in the orthorhombic phase, are illustrated in (c). 



the optimised, non-magnetic structures and is compa- 
rable to that obtained with the experimental structure. 
This is appropriate from a theoretical point of view as 
the phonon calculations should start from a fully relaxed 
structure. The optimised c-axis is within 0.5% of the 
experimental value and the As z-coordinate is close to 
the experimental value (see Table II). In addition, the 
calculated, orthorhombic distortion is strongest for the 
observed magnetic structure and closely matches the ex- 
perimental value. Magnetic ordering of spin states in Fe 
cations is therefore important for calculating FeAs bond- 
ing and giving the best, calculated vibrational density of 
states. The data in Table II also reveals that, from a 
structural point of view, the ferromagnetic case is simi- 
lar to the non-magnetic case. Also, including the correct 
magnetic ordering Mots^^ results in a volume expansion, 
which is of interest in light of the apparent pressure- 
induced superconductivity in the parent compounds of 
the FeAs-based superconductors^li^^. 



D. Effect of phonons on the electron density of 
states 



To gain insight into the effect of phonons on electronic 
structure, we have calculated the electronic DOS (eDOS) 
for frozen phonon patterns using the recently observed 
magnetic structure in the orthorhombic phase^i. The 
mean amplitude was chosen as ~ 0.05 A, which corre- 
sponds approximately to a single occupation of the mode. 
The results are shown in Figure 3(a) for the three Gamma 
point modes that most strongly perturb the equilibrium 
electronic structure. These modes are Raman active and 
are labelled B, C and E in Figure 2(c) and have frequen- 
cies 18, 26 and 32 meV. The main features of the eDOS 
for the unperturbed structure are in good agreement with 
other electronic structure calculations of the tetragonal 
phase^Si^i^. The bands close to the Fermi level and in 
particular the eDOS at the Fermi level depend strongly 
on the phonon distortion via the d-bands of Fe. As the 
Fermi level of the equilibrium structure is found at the 
minimum of a rather steep-walled valley of the eDOS, 
phonon modes that displace the Fermi level only slightly, 
create a large fluctuation of the electron density at the 
Fermi-level. 

In view of the importance of magnetic ordering on the 
equilibrium structure of BaFe2As2, we have also investi- 
gated the spin-resolved eDOS of the d-electrons of the 
Fe cation (see Figure 3(b)). Complex changes in the 
eDOS are captured in the magnetic moment of the Fe 
cation. The calculated magnetic moment for the equi- 
libiurm structure is 1.94 fis- This value is much larger 
than the observed magnetic moment ( 1 /is) as reported 
by other authors^^ and explained in terms of frustrat- 
ing magnetic interactions reducing the observed moment 
(see — and references therein) . On the computational 
side, standard settings, notably for band occupation at 
the Fermi level, tend to result in more localised and there- 



fore bigger moments. Settings could be tuned to give a re- 
duced magnetic moment but here we simply want to use 
the equilibrium magnetic moment as a reference value to 
quantify the effect of structural distortions. Mode B in- 
creases the moment slightly to 2.07 /is, mode C increases 
it to 2.55 iiB and mode E decreases it to 1.58 jib- From 
Figure 3(b) it is seen that modes B and C cause a decrease 
in the spin-down population and an increase in the spin- 
up population and therefore an increase in the magnetic 
moment. On the other hand, mode E causes a rocking 
of the spin population from the majority, up-channel to 
the minority down-channel, decreasing the magnetic mo- 
ment. Mode E tends to reduce the magnetic ordering. 
For all 3 modes, the Fermi level changes by less 0.5%. 



V. CONCLUSION AND PERSPECTIVES 

Small variations in the GDOS between the tetrago- 
nal and orthorhombic phases are observed in the new, 
high resolution, experimental data and these are gener- 
ally reproduced by the calculations. The variations in 
the distances between Fe atoms and between Fe and As 
atoms are less than 0.01 A^, revealing the sensitivity of 
phonons to subtle structural changes. Phonon calcula- 
tions performed on the experimental structure and on the 
fully relaxed structure, including the observed magnetic 
ordering, best match the experimental spectra. Non- 
magnetic and ferromagnetic, relaxed structures show big- 
ger discrepancies with the measured spectra. The key 
structural parameter is the z-coordinate of the As atom 
which governs the distortion of the Fe-As layer. Phonon 
modes that modulate the Fe-As layer cause the most sig- 
nificant changes in the eDOS at the Fermi level, which 
indicate the extent of electron- phonon coupling. In ad- 
dition the magnetic moment and spin-resolved eDOS of 
the Fe cation vary significantly with the phonon displace- 
ments, demonstrating spin-phonon coupling, which is rel- 
evant to the search for a composite mechanism for high 
temperature superconductivity that goes beyond simple, 
electron-phonon coupling. 

The main short-coming of the calculation is the inability 
to reproduce the observed, experimental peak at 20 meV. 
The calculation predicts in this region highly dispersed 
branches leading to a pseudo-gap in the density-of-states. 
Experimentally, the Q-dependence of this peak is similar 
to that of other phonon bands. There is, therefore, no 
experimental indication that we are dealing with a mag- 
netic excitation, especially at the temperatures of these 
measurements. In recent work by Reznik and coworkers^ 
combining ab initio, linear response, phonon calculations 
and inelastic X-ray scattering, the peak close at 20 mcV 
is observed experimentally but not reproduced computa- 
tionally. 

We have also tried to generate spectral intensity at 20 
meV within the framework of our phonon calculations. 
This is achieved by trying to fit inter-atomic force con- 
stants to change the frequencies of selected phonons at 
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FIG. 3: (Color online) Bottom panel (a): total eDOS of the unperturbed (Equilibrium), and the structurally perturbed 
orthorhombic phase of BaFe2As2. Total eDOS of spin-down channel is identical to the up-channel. Top panel (b): the 
corresponding spin-resolved eDOS of the d-electrons per Fe cation. Labels B, C and E indicate phonon-perturbed equilibrium 
structures (cf Figure 2(b) and Figure 2(c); B: 18 meV, C: 26 meV and E: 32 meV). 



particular wave-vectors. The only phonon that could be 
significantly modified in this way was a mode at 24 meV 
at the zone boundary (0,0,0.5), but the change in spectral 
intensity is small. This effect is however consistent with 
the sensitivity that we have demonstrated of the phonons 
to the z-coordinate of the As atom. The importance of 
the c-axis phonons in FeAs based superconductors has 
also been highlighted recently'^'^ . 

However, if we believe that modern DFT codes correctly 
reproduce the main spectral features in, structurally, 
rather simple systems like BaFe2As2, then we have no 
clear explanation for this peak within the phonon model 
implemented here. Therefore, any statement on the ori- 
gin of the 20 meV peak without further experiments 



has to remain speculative. Strong, non-conventional 
electron-phonon coupling has been proposed in the lit- 
erature^ for LaOFeAs. The possibility of a singularity 
close to 20 meV is evoked on the basis of gap openings at 
the Fermi surface. The same physical arguments should 
hold for BaFe2As2 as only the FeAs-layers are concerned. 
Finally we plan to extend our measurements and calcu- 
lations to other systems in this family of compounds, to 
determine whether the anomaly at 20 meV is a general 
feature or whether it is specific to the Ba system. In this 
context, we are motivated by initial work on the stron- 
tium system for which we obtain excellent agreement be- 
tween measured and calculated Raman frequencies. 
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